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Abstract 



Energy-dependent solar neutrino flux reduction caused by the Mikheyev- 
Smirnov-Wolfenstein (MSW) effect is applied to the Exact Parity Model. 
Several scenarios are possible, depending on the region of parameter space 
chosen. The interplay between intergenerational MSW transitions and vac- 
uum "intragenerational" ordinary-mirror neutrino oscillations is discussed. 
Expectations for the ratio of charged to neutral current event rates at the 
Sudbury Neutrino Observatory (SNO) are estimated. The implications of the 
various scenarios for the Boron neutrino energy spectrum and BOREXINO are 
briefly discussed. The consequences of MSW-induced solar neutrino depletion 
within the Exact Parity Model differ in interesting ways from the standard 
v e ^u,t an d v e <-> v s cases. The physical causes of these differences are 
determined. 



I. INTRODUCTION 

In the Exact Parity Model (EPM) fl[], parity is an exact symmetry of nature despite the 
V — A character of weak interactions. Exact parity symmetry is achieved by introducing 
parity or "mirror" partners for each of the standard model fermions, Higgs bosons and 
gauge bosons. In general, colour singlet and electromagnetically neutral particles in the 
standard sector mix with their corresponding mirror states, leading to possibly observable 
experimental effects. 
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One of the most interesting possibilities in this regard is mixing between ordinary and 
mirror neutrinos ||. In part, the EPM is an explicit theory featuring three effectively sterile 
light neutrino flavours in addition to z/ e , and v T . We shall denote the mirror neutrino 
flavours by is' e , v' and z/, where v' a is the parity partner of v a [a = e,fi,r). Exact parity 
invariance imposes a simple and nontrivial constraint on standard-mirror neutrino mixing: 
in the absence of intergenerational mixing, the mass eigenstate neutrinos must be maximal 
mixtures of ordinary and mirror neutrinos. This follows immediately from the requirement 
that parity eigenstates must also be eigenstates of the Hamiltonian when parity is an exact 
symmetry. The mass/parity eigenstates are given by 

Va± = — 71 — ' {) 

where v a ± — > ±(z/ Q= |-) c under a parity transformation. When intergenerational mixing is 
nonzero and CP violation absent, the mass eigenstates are simply linear combinations of 
the v a+ and, separately, the v a _: 

a 

^-=££/^a-, (2) 
a 

where i = 1, 2, 3 and U ia are unitary mixing matrices. Exact parity symmetry forbids 
mixing between positive and negative parity neutrinos in the vacuum.^ 

The neutrino sector of the EPM is of great interest because it can explain both the solar 
and atmospheric neutrino anomalies ||. The clearest case is provided by the atmospheric 
neutrino anomaly. (Note that we will consider the case of small intergenerational mixing 
in this paper, taking our cue from the almost diagonal Kobayashi-Maskawa matrix in the 
quark sector.) The observed anomalous value of the ratio R of //-like to e-like events strongly 
suggests that atmospheric muon neutrinos undergo large amplitude oscillations into another 
flavour. Large amplitude oscillations imply a large mixing angle between and another 
flavour. This is exactly what is provided for in the EPM through maximal <-> v' mixing. 
Furthermore, the anomalous zenith angle dependence for multi-GeV //-like events reported 
by SuperKamiokande provides strong independent evidence in favour of large amplitude 
oscillations of z/ M . The totality of atmospheric neutrino data is well explained by v' 
oscillations with Am\ +2 _ in the approximate range 

1(T 3 - Am 2 2+2 _/eV 2 - 10~ 2 , (3) 

where Am 2 , +2 _ is the squared mass difference between z/ 2+ and v 2 _ ||. Note that if we 
restrict the discussion to two-flavour oscillations, then the present data allows only two 



1 If the minimal standard model is extended by adding a mirror sector, then both neutrinos and 
mirror neutrinos are massless and unmixed. We do not consider this case because it is of little 
interest for neutrino phenomenology. 
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choices: the atmospheric neutrino problem is solved either by <-> v s oscillations (for 
which the EPM provides an explicit theory) or by <-> v T oscillations (see Refs. [||,|5| for 
a phenomenological study). In the future, these two alternatives may be experimentally 
distinguished through neutral current effects [§], upward through-going and stopping muon 
data |7||1, and long-baseline experiments. 

The solar neutrino problem also provides strong evidence in favour of the EPM neutrino 
sector. GALLEX || and SAGE [[LTJ observe a solar u e flux that is close to half of that 
expected from the standard solar model when neutrino oscillations are absent. A 50% u e 
flux reduction is exactly what is expected from the EPM due to maximal v e <->• v' e oscillations 
for the mass range 

10- 10 £ Am 2 1+1 _/eV 2 ~ 9 x 10" 4 , (4) 



where the upper limit is required for consistency with the CHOOZ bound [IT]. The other 
oscillation parameters are placed within the large region of parameter space where intergen- 
erational solar u e oscillations are unimportant. GALLEX and SAGE arguably provide the 
most unequivocal information regarding the nature of the solar neutrino problem. There 
are two reasons for this: First, theoretical calculations of the expected event rates are the 
most robust. Second, both detectors have been calibrated with respect to a neutrino source 
of known intensity. Kamiokande ||12|| , SuperKamiokande ]13|] and Homestake [14[] also pro- 



vide important information about solar neutrinos. All three of these experiments report a 
significant deficit of solar neutrinos, leading to a qualitatively consistent picture of solar v e 
depletion across the five experiments (see Table |). Furthermore, because of the different en- 
ergy thresholds of the experiments, a comparison of their results provides information about 
the energy- dependence of the solar neutrino flux depletion. Unfortunately, the precise sig- 
nificance of the information obtained from Kamiokande, SuperKamiokande and Homestake 
is less clear than for GALLEX and SAGE for the following reasons: (i) Predictions for the 
Boron neutrino flux vary significantly between different versions of the standard solar model, 
mainly because of an uncertain p + 7 Be — > 7 + 8 B cross-section. The precise value of the 
Boron neutrino flux deficit is therefore not as well established as one would wish, (ii) The 
pioneering Homestake experiment is still the only experiment that is especially sensitive to 
the mid-energy Beryllium neutrinos. Other experiments are needed in order to confirm their 
result. Fortunately, BOREXINO and the Iodine experiment will probe a similar part of the 
spectrum in the near future. They will either confirm or disconfirm the somewhat greater 
flux reduction reported by Homestake. 

The purpose of this paper is to study the range of possibilities for solar neutrino flux 
depletion provided for by the EPM, and to determine the implications of these possibilities 
for, in particular, the Sudbury Neutrino Observatory (SNO) [[TJJ. SNO will play a very 
important role in testing the EPM because of its ability to distinguish between solar v e «-» 
T oscillations and v e <-> v' a oscillations through its sensitivity to both charged and neutral 
current reactions. 

It is important to understand that the EPM supplies different solar neutrino outcomes 
in different regions of parameter space. Hitherto 0, work within the EPM has focussed 
on the simplest and most characteristic possibility: First, parameters are chosen so that 
intergenerational solar v e oscillations are unimportant. Maximal v e <-> v' oscillations in 
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the range of Eq.(^) then lead to an energy-independent 50% flux reduction compared to 
no-oscillation expectations. Furthermore, since v' e states are blind to the neutral current, 
this case leads to an expectation that SNO will measure the standard rate for charged 
current relative to neutral current events. This case is in many ways the most attractive 
possibility within the EPM, because it is extremely simple and because it most fully utilises 
the predictive power of the EPM: a 50% flux depletion is the direct result of maximal v e <-> v' e 
mixing which in turn is the direct result of exact parity invariance. 

However, this case does not reproduce the greater depletion of mid-energy neutrinos 
that is inferred from a comparison between the Homestake rate and the other measured 
rates. In this paper, we will explore regions of parameter space for the EPM that are 
different from that considered above and hitherto. There are two principal motivations for 
doing so. First, we want to identify those regions of parameter space that can provide a 
better fit to the totality of solar neutrino data than can the 50% flux reduction region. 
Second, since this study will necessarily involve solar v e oscillations into v' v T and v' T 
(as well as into z/), we will provide interesting predictions for the rate of charged current 
to neutral current events expected at SNO. Our study will essentially be an exploration of 
the Mikheyev-Smirnov-Wolfenstein (MSW) effect within the EPM. Note also that 

maximal v e <-> i/ oscillations can provide an energy-dependent flux reduction factor that 



fits all existing experiments well in the "just-so" regime fll9fl , 

Amj +1 _ ~5^8x W~ n eV 2 . (5) 

Thus an experimental confirmation of an active-sterile just-so scenario would also be con- 
sistent with the EPM for a tiny region of parameter space. 

Before commencing the analysis we should comment that Big Bang Nucleosynthesis 
poses a challenge for any model of light sterile neutrinos, due to the possible excitation of 
excess degrees of freedom during the relevant cosmological epoch. Fortunately, and indeed 
remarkably, sterile neutrino models can generally meet this challenge in full through the 
phenomenon of lepton asymmetry generation by active-sterile or active-mirror oscillations. 
For a complete discussion, see Ref. p0| . 

The rest of this paper is structured as follows: In Sec. |1| we provide a semi- quantitative 
analysis of the implications of the SuperKamiokande measurement for the ratio of charged 
to neutral current event rates at SNO. Section |T| deals with the mathematical formulation 
of the MSW effect on the EPM's underlying maximal mixing framework. We present various 
MSW solutions to the solar neutrino problem within in EPM in Sec. |V]. In Sec. |V|, we discuss 
the implications of the solutions for the Boron neutrino energy spectrum and the Beryllium 
flux, and reexamine the SNO charged to neutral current event rate, now constrained by five 
experiments. We demonstrate that consistency with the LSND result can be attained in 
Sec. 0. We conclude in Sec. |VTJ. 

II. FROM SUPERKAMIOKANDE TO SNO 

The Sudbury Neutrino Observatory (SNO) will determine whether solar i/ e 's oscillate 
into the other active flavours and v T or into sterile flavours or, to some level of sensitivity, 
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into a combination of active and sterile flavours. The last possibility is the generic prediction 
of the EPM, because mirror neutrinos are sterile with respect to ordinary weak interactions. 

In this section, we will outline the various solar neutrino outcomes possible in the EPM 
in different regions of parameter space. Our aim in this section is to estimate the ratio of 
charged to neutral current event rates that SNO will measure if the EPM is the correct 
theory of neutrino mixing. In this "warm-up" section, we will use as little theoretical input 
as possible in order to not obscure, by the technical complications of the MSW effect, the 
important phenomenological role played by the characteristic maximal u a <-» u' a oscillations 
of the EPM. Our minimal input will be: (i) that some type of MSW effect occurs (except 
for Case A below); (ii) that averaged maximal vacuum v a <-> 2/ oscillations occur when the 
oscillation length is much less than an astronomical unit; and, (iii) that SuperKamiokande 
has measured the correct depletion factor for Boron neutrinos. We will revisit this issue in 
Sec. [V], armed with more detailed information about the required v e survival probability. 
This will allow us to further constrain our estimation of the ratio of charged to neutrino 
current event rates expected at SNO. 

A. Case A: vacuum v e *-* v' e oscillations only 

Case A results from the parameter space region discussed in the Introduction and in 
previous papers ||. Parameters are chosen so that the only oscillation mode important for 
solar neutrinos is v e <-> u' e with the oscillation length set by Eq.(|J). If this case is correct, 
SNO will measure the standard value for the ratio of charged current to neutral current 
events. They should also, of course, confirm the substantial depletion of Boron neutrinos 
reported by Kamiokande and SuperKamiokande. 

As a subcase of Case A, another possibility is maximal v e «-> v' e oscillations in the "just- 
so" regime. The relevant parameters are given in Eq.fl5|). This energy-dependent subcase 
is essentially an active-sterile "just-so" scenario where the observed maximal mixing arises 
from exact parity invariance. In this picture, SNO will measure roughly the standard value 
for the ratio of charged to neutral current events. 

B. Cases employing the MSW effect 

Suppose intergenerational solar v e oscillations are now switched on by choosing a different 
point in parameter space. In order to obtain substantial intergenerational oscillations while 
simultaneously keeping the relevant mixing angles small, the MSW mechanism must be 
invoked. We will discuss the details of MSW transitions within the EPM in Sec. |T|. For 
the purposes of this section, we will merely suppose that MSW transitions exist and deplete 
solar neutrinos in an appropriate energy- dependent fashion. 

The interesting issue for SNO is the flavour content of the solar neutrino flux at the 
Earth. In general, MSW transitions will process some of the solar v e flux into second and 
third generation neutrinos and mirror neutrinos in the interior of the sun such that 

o (e, E) = (e, E) + (e', E) + (/x, E) + (/i', E) + (r, E) + (r', E), (6) 
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where (e, E) is the no-oscillation standard solar model flux of v e of energy E at the surface 
of the sun, while (p e (a,E) and 4> e (a',E) are the fluxes of v a and z/, respectively, at the 
surface of the sun. The equality in Eq. follows from flux conservation. The various fluxes 
on the right hand side of Eq. (^j) are given by 

^(a,E) = P°(E)<j ) °(e,E), (7) 

where P® a (E) is the matter- affected oscillation probability at the surface of the sun for 
v e *^ v a { a — oi,a/). Probability conservation requires that 

1 = P®m + P®*{E) + P°(E) + P e ® , (E) + Pf T {E) + P°(I2) (8) 

for each value of E. 

The crucial point can now be made: between the sun and the Earth, additional large 
amplitude vacuum oscillations will in general occur between the maximally mixed standard 
plus mirror pairs. Vacuum intergenerational oscillations will be small given our assumption 
of small intergenerational mixing. Provided that the oscillation lengths for v a u' a are 
smaller than an astronomical unit (i.e., where the corresponding squared mass difference is 
10 -10 eV 2 ), maximal vacuum oscillations will induce 



n a ,E)=n*',E) = * e{a ^ + 2 ^' E \ (9 ) 

where 0® denotes the flux at the Earth, and ( is a geometric factor due to the inverse square 
law. This is expressed in terms of oscillation probabilities as 

Pl{E) = Pg,(E) = P - (E) + P - /(jS) , (10) 

where the superscript © denotes oscillation probabilities at the Earth. The different possi- 
bilities for SNO now correspond to different oscillation length regimes for z/ a <-> v' a . 

Because strong evidence now exists for an atmospheric neutrino anomaly, we choose 
^ m 2+2- i n ^e range of Eq.(|3|). This means that Eqs.(0) and ([II]) certainly hold for a = \i. 



The different cases discussed below correspond to the four generic possibilities for Am 2 



and Am| +3 _ 



1. Case B 

Case B corresponds to the parameter choice 

Am 2 1+1 _, Am 2 3+3 _ Z 1(T 10 eV 2 , (11) 

so that Eqs.(||) and ( |10"D hold for a = e, r as well as for a = \i. Combining Eqs.@, (|9]) and 
(PI) we see that the total flux of active neutrinos at the Earth will be 

0> « ) = ^ = «|A (12) 
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that is, exactly half of the no-oscillation v e flux. This 50% flux reduction is a direct result 
of the maximal mixing constraint following from exact parity symmetry. This prediction 
implies that the SNO neutral current rate will be 50% of the no-oscillation expectation. 

To quantify expectations for SNO, we consider the rates for charged current and neutral 
current events given, respectively, by 



/■oo 

T cc = P®(E)<j)®(e,E)acc(E)dE, 
Je 

rod 

Tnc = / [P*{E) + P*(E) + P®(EM(e } E)a NC (E)dE } (13) 

JE 

where E is the energy threshold for SNO and <Tcc(E) 1<j nc (E)} is the charged [neutral] 
current cross-section. 



According to Eqs.(|), fllP]) , (p|) and (|13D , we see that the charged to neutral current rate 



divided by the no-oscillation expectation is given by 

Td = (r iv M = T — F ' ( 14) 

(1 CC/L NC)\0 J-CClo 

where the characteristic factor of two is just another way of expressing the 50% flux reduction 
of the sum of active flavours. SNO will measure T C c\osc, while Tc^lo depends on Boron 
neutrino flux predictions from the standard solar model. 

Equation (|T^) is an exact result. In order to obtain a precise prediction for r d , the 
energy-dependent survival probability must be known. However, a good estimation for 
can be obtained from the measured Boron neutrino flux at SuperKamiokande, because the 
energy threshold of SNO is similar to that of SuperKamiokande. The charged current event 
rate at SNO relative to the no-oscillation expectation should be approximately equal to the 
analogous quantity measured by SuperKamiokande. In order to use the SuperKamiokande 
measurement of flsx, where 

^ observed event rate . . 

"» S no-oscillation event rate ' (15) 

we have to correct for the small contribution that neutral current induced v^ )T e scattering 
makes to it. Using Eqs. (|T2D and ([14]) together with the relation a SK {v^ T e) ~ \o SK {y e e) 
between the relevant cross-sections at SuperKamiokande, we obtain 

r d ~ (16) 

for Case B. Using the information in Table |, and taking a 2a limit that incorporates both 
experimental and theoretical uncertainties, we get that 

tt SK ~ 0.25 -> 0.5. (17) 

The large range displayed here is mainly due to the significant theoretical uncertainty in the 
Boron neutrino flux. [Note also that we have focussed on SSM-BP (1995) only | 2l| . Other 



SSM calculations yield significantly different Boron neutrino fluxes f22fl .1 Using Eq. (|17j) , we 

get 
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r d ~0.4->l. (18) 

Note that for the upper extreme, where the Boron neutrino depletion is entirely due to 
averaged maximal v e <-> v' e oscillations, Case B becomes identical to Case A for Boron 
neutrinos. 

By way of comparison, the standard v e <-> v^ T expectation is 

r d ~ msK ~ 1 ~ 0.1 - 0.4. (19) 
5 

(Note that SuperKamiokande data only have been used to obtain this estimate. Tighter 
predictions are obtained when a survival probability consistent with all five solar neutrino 
experiments is used.) So, we see that the Case B range for covers all of the values between 
the ranges for the standard v e <-> v^ T and v e <-> v s solutions. A clear distinction between 
these three possibilities therefore seems to be provided for by r^, unless by bad luck the 
measured value turns out to be close to either 0.4 or I. 



2. Case C 

Case C is defined by the parameter choice 

Amj +1 _ Z 10- 11 eV 2 , Amj +3 _ Z lO" 10 eV 2 , (20) 

so that vacuum v e <-» v' e oscillations do not occur. In addition, both direct transitions of v e 
to u' e , and indirect transitions via the second and third generation flavours, are negligible 
within the sun. We can therefore set 

P°{E) = (21) 
to a very good level of approximation. In this case, 

r{active , E) = «h3±£i^l. (22) 

The total flux of active flavours is larger than for Case B given the absence of vacuum 
v e <-> v' e oscillations. 

The ratio of charged to neutral current rates for this case relative to no-oscillation rates 
is given by 

_o J^P^(E)^(e,E)a cc (E)dE f~J®(e, E)a NC (E)dE 
Td I^[l + P e ^EM(e,E)a N c(E)dEJ^J^e,E)a cc (E)dE^ 1 } 

where an explicit expression for the energy-dependent v e survival probability is required for 
an exact prediction. In principle, this has to be done on a case by case basis. 

An approximate indication of the likely outcomes is obtained by neglecting the energy- 
dependence to obtain 
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2(P®) 

rd 1 + <P®>' 1 j 

where (. . .) denotes an average. Taking the SuperKamiokande measurement of Qsk, and 
correcting for neutral current effects using Eq. (|22|) , we obtain 



(PS) ~ H5£pi. (25) 

For in the range of Eq.([P7D, this implies that 

r d ~ 0.3 -> 0.6. (26) 

This case should be clearly distinguishable from Case A. It is not distinguishable from Case 
B or from the standard MSW v e — > v^ T scenario on the basis of r d alone. We discuss this 
issue further in Sees. M and [VII. 



3. Case D 

Case D corresponds to 

Am 2 1+1 _ £ 10- 10 eV\ Am 2 3+3 _ & lO' 11 eV 2 . (27) 

No especially interesting predictions can be made in this case without further information. 
For instance, if the MSW partners of u e are and v' then this case reduces to Case B. If, 
on the other hand, the MSW partners of v e are v T and u' T , then, in the energy-independent 
approximation, 

Ti ~ {pg> +>«?)' (28) 

For (P®) = 0.25, probability conservation at the surface of the sun requires < (P® ) < 0.5, 
leading to r d > 0.33. The lower bound on r d increases with (P®), going to 1 as (P®) 
approaches 0.5. 



4- Case E 

Finally, Case E corresponds to 

Am 2 1+1 _, Am 2 3+3 _ ~ 10" n e1/ 2 . (29) 

Again, more information is needed in this case in order to obtain predictions. If the MSW 
partners of v e are and v' then this case reduces to Case C (and is also similar to the 
scheme analysed in Ref. 0]). If, on the other hand, the MSW partners of v e are v T and 
z/, then this case is intermediate between the standard v e ^ v^ T and u e — > v s scenarios 
because vacuum oscillations play a negligible role. 
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III. MATHEMATICAL FORMULATION OF THE MSW EFFECT IN THE EPM 



The EPM predicts an energy-independent 50% reduction of the solar neutrino flux by 
maximal vacuum v e «-> u' e oscillations that is in good quantitative agreement with experi- 
ments primarily sensitive to low-energy and to high-energy neutrinos. The significantly lower 
event rate measured by Homestake, however, calls for further suppression of the mid-energy 
flux. Preferential energy-dependent depletion can be achieved via the MSW mechanism by 
restricting the relevant intergenerational squared mass difference to |P7| 



10~ 8 ~ Am 2 /eV^ - 



t2 < 



10" 



and jHJ, 



• 2 o > 
sin 2r] ~ 



10" 



(30) 



(31) 



where r\ parameterises the mixing between the corresponding neutrino states that take part 
in resonant conversion. 

In the following analysis, a standard neutrino mass hierarchy, i.e., 



2 . 2 ^ 2 
m 3± > m 2± > m 1± , 



(32) 



is assumed such that the MSW partners of v e are and v' . Under the assumption of small 
intergenerational mixing, the contributions of v T and v' T towards the v e survival probability at 
the Earth through vacuum oscillations is negligibly small. It then suffices to consider only 
the interactions between the first two generations, though the forthcoming mathematical 
treatments can be easily generalised to include the third generation. The remaining two- 
generation system thus consists of four neutrino states, where the transformation between 
the weak and mass/parity eigenstates is given by Eqs.([TJ) and (Q). Explicitly, 



(Pe\ 




( 




C e 


S<f> 


Se\ 




(K- 
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1 




—C^ 
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S e 
















—Sg 
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C e 
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—C^ 


C e ) 




\ "2+ 


) 



(33) 



that are responsible 



where 9 and parameterise the two 2x2 unitary matrices Uf a in Eq. 
for the respective mixing of positive and negative parity eigenstates, and — | < 6 
Exact parity symmetry thereby reduces a nominally six-angle problem (if CP is conserved) to 
a two-angle task. (A generic 4x4 orthogonal matrix consists of six independent parameters.) 

Note that, at this stage, we do not make any assumptions regarding the signs of Am^ +1 _ 
and Am2 +2 _, and Eq.(|33|) does not imply in any way that v\ + (^2+) is heavier than z^i_ (^2-)- 
This is because, firstly, we have no prior reasons for doing so. Secondly, since neutrino states 
of unlike parity do not mix in a gauge theoretic sense (i.e., the Lagrangian of the EPM in 
vacuum does not contain parity- violating terms such as mz/ 1+ z/ 2 _ for z/ 1+ <-> v 2 - @)> we 
would expect the effective mixing of like-parity (such as u 1+ <-> u 2+ ) and of unlike-parity 
eigenstates to experience different forms of matter enhancement. Ours being a two-angle 
problem renders the quantification of this difference a relatively simple task. Henceforth, 
we shall denote the heavier (lighter) of z/ 1± and of u 2 ± as v lh and u 2h {y u and v 2 i) such that 
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(34) 



in order to keep the analysis as general as possible. 

The problem now becomes one of solving the Schrodinger equation 
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(35) 
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ml e and 



where U is the mixing matrix in Eq.(^), E is the neutrino energy, and m 
m 2h are the squared masses of the mass/parity eigenstates uu, v lh , v 2 t and v 2 h respectively 
The interaction terms for v P and v.. 



[v' e and v' are inert) are 



A„ 



Acc + Anc, 

AnCi 



(36) 



where CC stands for charged current, NC for neutral current and 

A cc = 2V2G F EN e (x), 
A NC = -V2G F EN n (x), 



(37) 



where Gf is the Fermi constant, N e (x) the electron number density at position x in the 
neutrino's path, and N n (x) the neutron number density. If U m is a density-dependent 
unitary transformation that puts the total Hamiltonian 7i in an instantaneous mass basis 
vT" such that 



*a = £ U*v?, a = e, e', ft, //, % = U, lh, 2£, 2h, 



(38) 



the Schrodinger equation in Eq.([35|) can be rewritten as: 



dx 



11, 



v 

V 2i 
\^2H 



re 



V, 



'lh 

,771 
2/' 



1 

2E 



( 



u m -\n - i^-)u r > 

dx 









m 



ih 



(x) 











m 2 2e {x) 




1/1 

,771 

2( 



V < h J 









-u 



m—1 



dx 



m 



2i, 



V 

V. 



It 
rn 
lh 

m 
2f 
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where m 2 £ { 



x 



m lh( x ), m li( x ) and m lh{ x , 
neous mass eigenstates z/^, z^, and respectively. Given the initial conditions 



are the squared mass eigenvalues of the instanta- 



,(Xi) = 1, v' e {xi) = V^Xi) = v'Jxi) = 0, 



(40) 



where Xi is the u e production position, the probability that a v e produced in the sun will be 
detected on Earth is expressed as 



^C/ e 7(^)exp 



1,3 



-i l Xf H%dx 



(41) 



where i, j = l£, lh, 21, 2h, the exponential of the integral from X{ to the detection position 
Xf is the solution to Eq. (|59"D , U m (xi) is the density-dependent mixing matrix at the point 
of v e production and we have chosen U real (assuming that CP is conserved). The term 
exp — i f^f TQjdx represents the amplitude of a transition from v™ to uj 1 and vice versa. 

For future reference, note that for a 2v system, after phase-averaging, the u e survival 
probability is given by [E4j] 



P®{E)\ 2v = - + (^ - P R ) cos2r) m (xi) cos2r/ = P 2u msw, 



(42) 



where rj is the vacuum mixing angle, rj m ( x i) the matter mixing angle at the u e production 
position, and Pr the level-crossing probability evaluated at resonance. [A common practice 
is to multiply the Pr term by a step function, 6(E — Ea), where Ea is the minimum energy a 
neutrino produced at X{ must possess for a resonance to take place inside the sun p5f . It shall 
be omitted for the convenience of typesetting.] We shall refer to this 2v survival probability 
as the MSW transition probability PivMsw- Bearing in mind that the same expression can 
be obtained by considering probabilities instead of amplitudes, we shall adopt the same 
classical attitude for the rest of the analysis. 

An analysis involving two effectively sterile neutrinos approximately maximally mixed 
with v e and respectively was carried out in Ref. fZEfl . The "intragenerational" mass 
differences were assumed to be much smaller than the intergenerational mass difference, the 
latter of which was responsible for a v e <-» MSW resonance. Our case differs in that the 
explanation of the atmospheric neutrino anomaly by maximal vacuum <-> u' oscillations 
requires the squared mass difference between v 2 + and z/ 2 _ to lie within the range of Eq.(^). 
These masses are much larger than the MSW masses in Eq. (|30D , i.e., 



Ami, 



_ 2 _ > Am 2 21 , 



(43) 



where Am 21 = m 2£ — m\ h . On the other hand, the squared mass difference between 
Vig and vn t is constrained only by an experimental upper bound of 9 x lO'^eV 2 from 
Eq.(|]). Thus, analytically a priori well-approximated MSW solutions exist for three dis- 
tinct neutrino mass hierarchies, which we shall denote as Cases Bl, B2, and C respectively: 
Case Bl: Am\ +l _, Aml +2 _ > Am| 1; 

Case B2: Am 2 2+2 _ > Am 2 2l > Am 2 +l _ ~ 10- 10 eV 2 , 

> Am 2 2l > Am 2 1+1 _, Am 2 1+1 _ ~ lQ- u eV 2 . 



Case C: Am 2+2 _ 
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These somewhat playful labels are chosen for consistency with Sec. H, i.e., the generic pre- 
dictions for SNO for Case B in Sec. H hold for both neutrino mass hierarchies defined in 
Cases Bl and B2, and similarly for Case C. In the following subsections, we shall derive the 
u e survival probability for each case. 

A. Case Bl 

Case Bl assumes the following neutrino mass hierarchy: 



AmL_,Am2 +2 __ ^> Am 



21 • 



In this analysis, we make one further assumption that 



Am 



1+1- 



l(T 3 eV 2 



(44) 



(45) 



for simplicity. The squared masses of the instantaneous mass eigenstates for the relevant 
solar densities are shown in Fig. [y. 

We identify the point R, at which v™ h and vT^ almost cross, as an intergenerational MSW 
resonance. Large Am\ +1 _ and Am\ +2 _ (compared with Am^) ensure that u e <-> v' e and 
Ufj, <-> u f oscillations remain close to maximal in its vicinity. Consequently, matter effects 
are most strongly felt by vxh arid u 2 e, leading to the instantaneous mass eigenstates v™ h and 
v%l bearing little resemblance to their vacuum counterparts. The evolution of vu and v 2 h 
near R, on the other hand, is only slightly affected by matter, so that 



v 



V. 



u 

III 
2h 



(46) 



Provided Eq. (flop is satisfied, Eq.(fP]) will continue to hold at densities p ^> pn in the sun. 
Thus, to a very good level of approximation, 



u. 



e2h 



U e 2h, 



(47) 



which are virtually density- independent. [Note that if Amf +1 _ satisfies Eq.(f4"4]) but is, at 
the same time, sufficiently small (say, ~ 10~ 5 eV 2 ), matter effects can cause v e and u' e to 
depart from their mutual maximal mixing at p ^> pr. U'^^Xi) becomes density-dependent 
and may be drastically different from its vacuum counterpart. The quantification of this 
effect is relatively simple. For our purposes, however, we shall not consider it here.] Thus, 



'u 



uu) and v. 



2h 



i>2h) decouple from the system and evolve adiabatically such that 



after phase-averaging, we may write, following the procedures in Ref. fl27| , the v e survival 
probability as 



P^E) = \U eU \ 4 + \U e2h \ 4 + 



lUeiel 4 + \U e 2h\ 



i,j=lh,2£ L Jxi 

(1 — |^el£| 2 — \U e 2h\ 2 ) 2 P2vMSW, 



u, 



(■J 



(48) 
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where 



2vMSW 



E 

i,j=lh,2i 



TTm2 I Jjm2 
U elh "T u - 



exp 



1 e2l 



u, 



'-,) 



c2f 



(49) 



and U%1 + U, 



Tm2 

' e2t 



u 2 elh + u? 



2( 



by unitarity. The physical interpretation of Eqs. (ESI) and 



U. 



eU\ 



and 



follows simply. At well above the resonance density, a fraction \U^ e (xi)\ 2 

another \U™ h (xi)\ 2 (~ \U e 2h\ 2 ) of the i/ e 's produced at Xj populate the instantaneous mass 

1 



' e2h 

eigenstates v^} 

\Ue2H?~l 

where 



v\i ~ |) and v™ h (~ v 2 h ~ 0) respectively. The remaining (1 — \\J e \i\ — 
is distributed in v™ h and v™ t in a density- dependent ratio, cos 2 1r\ m : sin 2 2?] m , 



cos rj Ti 



jjm 
u elh 



Tjm2 I ffm2 



sin 7] ri 



jjm 

u e2l 



Tjm2 I Tjm2 

u elh u e2l 



(50) 



and 7] m ~ | for p ^> p^. So the z/ e 's produced in the region p ^> p/? near the centre of the 
sun populate u^}, u^l, vJJ} and v™ h in the approximate ratio | : 0| : 0. The half residing in 
v^ g participate in resonant conversion at R, while the other half in u^} (~ vu) propagate 
adiabatically to the surface of the sun without passing through a resonance. 

To study the intergenerational MSW resonance, we may treat the Vxh v 2i subsystem 
as forming two orthogonal pseudo-weak eigenstates, v a and h>b, that convert resonantly into 
each other at R, i.e., 







f u u \ 




(I 












( Vu \ 




= n( v ) 


V\h 








COST) 


sin rj 












V21 




— sin 77 


COST] 




V21 


\v 2h J 




V v 2h J 




\o 








I J 




V V2h J 



(51) 



where 



cos T] 



U, 



elh 



TJ 2 4- TJ 2 

u elh + u e2l 



sin 77 



u, 



<2I 



u elh + u e2l 



(52) 



are the vacuum counterparts of the parameters in Eq. (|50|) . With this parameterisation, we 
may rewrite the mixing matrix U as 



U = T (other s)TZ(rj), 



(53) 



where T(others) is a unitary matrix responsible for other mixing modes. Since vih «-> vu 
is the only matter-enhanced mixing mode, with the matter mixing angle r] m taking on the 
maximal value | at resonance, we may approximate U m as 



U m = T (other s)K(r) n 



(54) 



in the same manner that is adopted in the study of 3v schemes [ 25], 2"B|| . Equation ( [54] ) then 
allows us to recast the integrand TV£j in Eq.(|49|) into 
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rr = u m - l (rt-i^)u m 

dx 

= ll- 1 (r 1rn )(H e f f -i—)K(r 1m ), (55) 



where 



Heff = n(7 1 )u- 1 nun- 1 (7 1 ) 

= ^[n( v )Mn-\ v ) + n^u-'n^un- 1 ^)}, (56) 

by Eqs.(|35|) and (|39|). The 2x2 submatrix TC e ffij, where i, j = lh, 2£ (i.e., the 2 — 3 sector 
of Tl e ff), is thus the effective Hamiltonian that governs the evolution of the pseudo-weak 
eigenstates u a and v^. Subtracting a common phase, TC e ffij takes on the form 

J_ / -Am^ cos2?7 + A eff Am| sin 2r) + A ind \ , . 

where A e ff, the effective density, is 

A ef f = A cc (U 2 elh + U 2 e2l ) 

+ A NC [coa 2r}(U 2 lh - U 2 2e + U 2 lh - U 2 2e ) + 2 sin 2r](U elh U e2£ + t/^C/^)], (58) 

and 

A ind = A NC [2 cos 2r}(U elh U e2e + U^U^e) ~ sin 2??([/ e 2 1/l - U 2 2e + E/J lh - U 2 2i )], (59) 

where the subscript ind stands for induced. The physical significance of this term will be 
discussed in due course. The quantities A e ff and A in d evaluated for various combinations of 
Uih and v 2 i in the EPM are shown in Table [TT]. The effective Hamiltonian Tl e ffij is analogous 
to that for a standard 2v system. The solution to 



d 



will thus give us the term P^vmsw that is equivalent to the MSW transition probability given 
by Eq. fl42|) , with a squared mass difference of Am^ and mixing angle 77 defined in Eq. fl52|) 
in a medium of effective density A e ff given by Eq.(|58"l) (plus some non-standard features to 
be discussed). 

We now look at the 2v subsystem more closely. In general, the effective density A e ff 
in Eq.([58|) contains both charged and neutral current interaction terms, though the latter's 
contribution is negligible if intergenerational mixing is small, i.e., 

Aeff « \a cc , (61) 

for \U e ih |, \U n2 i\ ~ ^= and \U e2 (\, \U^ih\ ~ 0. Equation (|61|), in turn, supports an approxi- 
mate resonance condition given by 
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.4 



cc 



2Am 2 21 



cos 2^, 



(62) 



which may be rearranged into a more illuminating form: 

Amlx cos 2r] 



V2G F N e 



Xi I 



(63) 



where the subscript A stands for adiabatic and N e (xi) is the electron number density at the 
v e production position. The quantity determines the location of the adiabatic edge of 
the MSW transition probability in the limit of small intergenerational mixing, such that for 
all neutrinos produced at Xi, only the ones with energy E > will be resonantly converted 
into other species. Comparing Eq. (|63"D to its counterpart in the standard u e <-» T scenario 
(where A e ff = Acc), our naturally smaller effective density automatically puts the adiabatic 
edge at twice the energy of the latter for a given Am^ and rj. Equations (|6"ID, ( [52]) and ( |63| ) 
are exact and independent of intergenerational mixing (provided it is nonzero), according to 
Table [TT] if (i) v\k and v 2 i are both positive or both negative parity eigenstates, or (ii) 9 = (ft 
for all possible combinations of v\h and v 2 i- Nonetheless, while we are not considering large 9 
and cases, Eq. fl6"3|) will locate the adiabatic edge with sufficient accuracy regardless of the 
exact identities of v\h and v 2 t for the present analysis in the limit of small intergenerational 
mixing. Taking \U e ig\ ~ 775 and \U e 2h\ ~ 0, the v e survival probability for Case Bl in Eq. 
is well approximated by 



Pf e {E) 



(64) 



Equation (|64[ ) is plotted in Fig. ||, juxtaposed with the respective survival probabilities for 

cases evaluated for the same oscillation parameters for 



and v e 



v. 



the standard v e 
comparison. 

The second non-standard feature is the presence of density-dependent terms, A in d, in 
the off-diagonal elements of the effective Hamiltonian in Eq.(^), representing some form 
of matter-induced mixing similar to that discussed in Ref. 0]. This matter-induced mix- 
ing manifests itself primarily in the non-adiabatic high-energy end of the MSW transition 
probability. Explicitly, if we write the effective Hamiltonian as 



-£(x) k(x) 
k{x) £{x) 



n *ff = \ XV 'c)Z\ ) ' ( 65 ) 
assuming a linear density profile, the level-crossing probability may be written as p9|j30[ 

P R = exp 

where the adiabaticity parameter j R is 



TT 



"2 7 « 



(66) 



IR 



, k 2 {x) 



dt(x) 
dx 



(67) 



X=X R 



evaluated at resonance; Pr's dependence on the matter-induced mixing term A ind is obvious. 
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In the context of the EPM, the extent to which matter-induced mixing affects the non- 
adiabatic conversion of solar neutrinos depends largely on the identities of V\h and z/ 2 ^. In 
particular, the mixing of like-parity and of unlike-parity eigenstates receive considerably 
different forms of matter enhancement. With reference to Table m if V\h and v^t are like- 
parity eigenstates, A in d vanishes exactly, leaving behind in Eq.(|57|) the standard vacuum 
parameters Am^ and rj, where rj is now replaced with 9 or <fi for v\ + <-» z/ 2+ and v\- <-> z/ 2 _ 
respectively. Thus, parity-conserving, direct mixing modes are enhanced naturally by matter 
effects in a familiar resonant fashion. In addition, the resonant enhancement of one mixing 
mode is completely independent of the other, that is, if the resonant mode is v\ + <-» z/ 2+ 
where 9 is the mixing angle responsible, <j) does not enter the scene. 

In the EPM, the apparent mixing of unlike-parity eigenstates in vacuum is an observa- 
tional effect due to mixing through other parity-conserving channels. In matter, apparent 
parity-violating mixing is, to some extent, conjured up by matter, as suggested by the gen- 
eral presence of an A ini i term for the matter-enhanced effective mixing of V\_ «-> z/ 2+ and 
u% + «-> v 2 - respectively. Furthermore, the strength of this matter-induced mixing is depen- 
dent on the relative amplitude of the 9 and (f) modes. An inspection of Table [TI| reveals that, 
depending on the sign of sin(# — </>), matter-induced mixing may enhance or suppress the 
non-adiabatic conversion of u e by decreasing or increasing respectively the matter oscillation 
length at resonance. The magnitude of this matter-induced mixing is, in part, controlled 
by the neutron density at resonance but is most severe when 9 and (f) differ significantly. 
Consider the case of effective v\- «-> z/ 2+ mixing. The contribution from matter-induced 
mixing relative to vacuum mixing is represented by the ratio 



A 



A NC cos 2r\ sin( 



Am 21 sin 2r\ 



X=X R 



(68) 



evaluated at resonance. Using the resonance condition in Eq . (p2[) , we find that 



A 



N n (x R ) 
N e (x R ) 



cos 2 2rj sin(# — 



sin 2r] 



1 N n (x R ) 

2 NJx R ) 



1 - 



tan< 



tan# 



(69) 



by various relations in Table || to first order in sin?]. Matter-induced and vacuum mixing 
are comparable if A w 1. Given that the electron number density is some two to six times 
the neutron number density in the interior of the sun [21], this corresponds to 



t&n6 NJxr) 

— r < 2 v ; + 1 

tan6> N n (x R ) 
if 9 and are in the same quadrant, or 



5 -> 13, 



(70) 



tan0 
tan^ 



< 1 -2 



N e {x R ) 
N n (x R ) 



-11 -> -3, 



(71) 



if 9 and <fi are in different quadrants, in order for matter-induced mixing to be recessive. 
The most extreme scenario is when the 9 mode is completely absent, such that sin 77 = 
and matter-induced mixing completely dominates. The consequential shift of the adiabatic 
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edge in the MSW transition probability is negligible. On the non-adiabatic side, assuming 
a linear density profile, the level-crossing probability is determined by 



R 



cxp 



exp 



4 E 



1 dA eff 



dx 



71 Amji sin 2 ! 



X=X R 

I Awe I 



4 E 


1 dA eff 




A ef f dx 



2Aml l cos 



Comparing this with its counterpart in the case where 9 
exactly, see Table [TT]), that is, 



P R = exp 



7i Am^i sin 2 2ip 



4 E cos 2ip 



1 dA eff 



dx 



(72) 

x=x R 

ip (such that A in d vanishes 

(73) 



X=X R 



we observe that approximate agreement between Eqs. (|72|) and ([73]) in the small intergener- 
ational mixing limit requires (for the same Am^), 



(2A 



rn 



sm 



21 x sin 2 2^ 



'2N, 



\ I Awe I / 
(16 -> 144) sin 2 293. 



N n (x R ) J 



sin 2yj 



(74) 



Thus, if we fit the z/ e survival probability in Eq.(p4|) to experimental data for the 9 = case 
(so that the identities of v lh and v 2 e do not matter), we know automatically from Eqs.flBTD 
to ( |63| ) that approximately the same Am^ will account for v\_ <-> z/ 2 + with = 0, while the 
mixing required, sin 2 20, is some 16 to 144 times that for 9 = <fi, under the assumption of 
small intergenerational mixing according to Eq.([74]). 

On the other hand, if is set to zero and 9 allowed to vary for v\- <-> U2+, we see 
from Table [TJ that both vacuum and matter-induced mixing contribute to the level-crossing 
probability, i.e., 



R 



cxp 



exp 



71 



Am 



21 



(sin 2t] + 



4^h cos 2ri sin 9 

Ami, 1 



4 E 



71 



cos2?7 
Am 2 x sin 2 29 ( x 



1 dA eff 



^4 cos 2 



2- 



<±t: 



X=X R 



4 £ 



1 rfA e/J 



dx 



cos 4 (9(1 + sin 2 



(75) 



where we have used various relations in Table [TJ. Hence, given 9 and 0's minute effects 
on the adiabatic edge and thus the fitted Am^, if sin 2 2 95 fits the data for the case where 
9 = 0, the case «-> z/ 2 + with = will be well described by the mixing parameter sin 2 29, 
which is approximately related to sin 2 2<p, assuming small intergenerational mixing, in the 
following manner: 
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sin 2 26 



1 - 



\A N c\ 
2AmZ, 



sin 2 2ip 



= (1.2 -> 1.8) sin 2 2y>. 
The same analysis applies to the case of v\ + <->• z^- • 



Nn(x R ) 

2N e {x R ) . 



sin 



(76) 



B. Case B2 



Case B2 corresponds to 



Am 2 +2 _ > Am 2 x > Am| +1 _ ~ 10~ w eV 2 . 



(77) 



The squared masses of the instantaneous mass eigenstates for the relevant solar densities 
are shown in Fig. [3[ As in Case Bl, the evolution of u%l is virtually density- independent 
and thus adiabatic due to a large Am\ +2 _. To a very good approximation, 



and 



'2h 



u e2h 



"2h, 



u, 



e2h- 



(78) 



(79) 



The decoupling of v 2 h renders the v e survival probability into the form: 



P®(E) = \U e2h \* + 



i,j=U,lh,2£ 
\U e 2h\ 4 + (1 — \U e 2h\ 2 ) 2 P , ivMSW-, 



u, 



where 

P?,vMSW 



E 



U m 



~iK2t x/U£ + U& + U% 



exp 



rx f 

-1 I HTjdx 



U, 



r/2 1 TJ2 , 112 
U ell ^ U elh T U e2E 



(80) 



(81) 



and U™ 2 + U™ 2 h + Ur 2 



eU ' ^elh ~ ^ e2l ~ w ell 

u 2h ) state as containing a fraction \U I 



eih + ^e2e unitarity. If we regard the decoupled v\\ 



m 12 

e2h\ 



\U e 2h\ ~ 0) of the original v e population, 
the other (1 — |t/ e 2?t| 2 ~ 1) is thus distributed in the remaining three states. The splittings 
between these states lie, by assumption, within the MSW range given in Eq.(^), forming 
a Zv subsystem which undergoes, technically, two resonances Rh and Rl- In this manner, 
P-ivMsw is equivalent to the v e survival probability for a standard v e <-> <-» v T system 
with m 2 ^> m 2 ~ m 2 P5| , P8| , |5HJ52| (plus some non-standard features due to the presence of 
sterile neutrinos). Contrary to the standard 3u system where the spotlight is on Rl, and 
Rh occurs at too high a density to be relevant, our focus is on Rh- (We do not consider 
Rl because the close encounter of the mass eigenvalues of i/JJ and v™ h in Fig. [| represents 
maximal conversion of v e into v' e . This happens in vacuum, necessarily adiabatically.) 
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Standard 3u systems have been studied extensively [fjl|,^,[31],[3^] . Following from these 
analyses, we parameterise the following mixing angles: 



COS^n 



U, 



eli 



sin ip v 



jjm 
u elh 



'elh 



COS T] ri 



TTm2 I TTm2 
eli ' elh 
m2 T TTm2 T TTm2 ' SmT lm 
ell i* ^elh i* u e2l 



m2 _j_ Tjm2 
' u elh 

U e2l 



Tjm2 
U elt 



Tjm2 , rr 



m2 

r2t 



(82) 



where the subscripts and superscripts m denote their density-dependent counterparts. The 
angle ip m describes (approximately) the mixing of v e and z/, and takes on a value of ip m = 
ip = j in vacuum. This mixing mode is strongly suppressed at high densities because of the 
small vacuum splitting between v u and V\k (so that ip m — > | as p — > oo). Consequently, the 
z/i state that takes part in the intergenerational MSW resonance Rh at density p ~ Am^ 
resembles neither z/^ nor z^ but is, asymptotically, some approximately maximal linear 
combinations of the two states, which we denote as v\ x . Its orthogonal state, Vx y , is thus 
the asymptotic form of for the relevant solar densities, that is, 



V2t 
V V2h ) 



( VU \ 

V21 

V "2h J 



I cos ip sin ij) \ 

— sin if> cos ip 

10 

V 1/ 



( Vu \ 

V21 

\ V2h ) 



(83) 



Note that for small intergenerational mixing, v lx ~ v e and i>\ y ~ v e by Eqs. (p2[) and (p3| ). 
At density p pr h , we may treat the 4z/ system as consisting of two parity eigenstates 
V21 and u 2 h, and two pseudo-mass eigenstates V\ x and z/ig, of indefinite parity. 

The angle 77 describes the apparent mixing of v\ x and v^t in vacuum, which is minimal as 
inferred from Eqs. (|3"3"|) and ( |52"D assuming small intergenerational mixing. However, in the 
proximity of Rh (at p ~ Amf-J, while u ly and z/ 2 /i propagate virtually density- independently, 
matter effects rotate r\ m from its vacuum value 77 through | to ~, and thereby modify the 
evolution of V\ x and v<it dramatically. If we regard V\ x and z/ 2 £ as forming two orthogonal 
pseudo-weak eigenstates u a and v h that convert resonantly into each other at Rh, i.e., 
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cos 77 sin ip sin rj \ 
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V 



sin 77 sin ip cost] 
1/ 



Vlh 

V21 

V V2h ) 



34) 



we may follow a similar procedure to Eqs.(^) and (p4j) and approximate [7 m as 

tT" = t/0- 1 (^)5- 1 (r/) l S(7 /m )0(^ m ), 



(85) 



such that the effective Hamiltonian that governs the evolution of v a and Uf, in the vicinity 
of Rh is given by the 1 — 3 sector of 
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n eff = ^[s^o^mo-^s- 1 ^) + s^oim-'nintuo- 1 ^- 1 ^)} 

i-zseetor \ ( -Am 2 eff cos 2r) + A eff Am 2 eff sin 2 V + A ind \ cQnstant m 
4\ Am 2 eff sm2r] + A ind Am 2 eff cos 2rj - A eff J 



where 



Am^j = ml e — -{m 2 e + m\ h — AmJ +1 _ cos2^) = m 2 ^ — m\. (87) 



The term Am^j represents the effective vacuum squared mass difference responsible for 
the resonance, and m\ is the averaged squared mass of the V\ states, m\ = ^(m 2 £ + m"l h ). 
The magnitudes of Am 2 ^ and Am^, where Am^ = m 2l — m 2 h , are virtually identical. 
Henceforth, we shall replace Am 2 ^ with Am^ whenever the former is encountered for 
convenience. The quantities A e ff and A ind are given by 



A efS = A cc {U 2 eU + U 2 lh + U 2 2e ) 



+ A nc {cos2t][(cos i/;U e u + sm^U elh ) 2 - U 2 2£ + (cos^U^u + simpU^h) 2 - U 2 
+ 2sm2r}[U e 2e(cosil)U e u + smipU elh ) + U^eicos ipU^u + sin ipU^h)]}, 

and 

A^d = A N c{2cos2r][U e2e (cos^U eU + sinipU elh ) + U^coaipUpu + sin^f/^)] 

- sin2?7[(cosV'f/ e i£ + smtpU elh ) 2 - U 2 2£ + (cos^£/ M i£ + sinipU^h) 2 - U 2 2£ ]}, (89) 

respectively. Table fT| shows A e ff and A ind evaluated for various combinations of uu, V\u 
and V2t- 



Taking the well-established 3z/ survival probability from Ref. [25] and setting the level- 
crossing probability at Rl to zero, we may immediately write down P^msw as 

PzvMsw = sin 2 f] m {xi) sin 2 r] + cos 2 r) m (xi) cos 2 7?[sm 2 if) m (xi) sin 2 %j) + cos 2 if> m (xi) cos 2 ip] 

- P RH [sin 2 r] m (xi) - cos 2 rj m (xi) sin 2 Vv(^i)] (sin 2 77 - cos 2 rj sin 2 ip) , (90) 

where Pr h is the level-crossing probability at resonance Rh calculated from TC e ff in Eq.(^). 
Furthermore, because of the strong suppression of u e <-> ^ oscillations in most parts of the 
sun, the corresponding matter mixing angle ip m is close to |, thereby reducing Eq.(^) to 

P-ivKisw = sin 2 r] m (xi) sin 2 r? + cos 2 T] m (xi) cos 2 r? sin 2 *0 

- P RH [sin 2 ?7 m (xi) - cos 2 ?7 m (a;i)](sin 2 ?7 - cos 2 77 sin 2 t/>). (91) 

Equation (|9l|) does not present itself in the most illuminating for the purpose of comparison. 
In the first instance, it does not, superficially, resemble the familiar expression for the stan- 
dard 2v MSW transition probability Pi v msw m Eq.(^). However, putting it into context, 
the overall v e survival probability P®(E) for Case B2 in Eq. (|80|) can be recast into 

PZ{E) = \U e2h \ A + (1 - \U e2h \ 2 )(l - \U eU \ 2 - \U e2h \ 2 )P 2uM sw, (92) 

where 
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P&MSW = ^ + (9 " P Rh) COS 2 Vm{Xi) cos 2w, (93) 



and 



COSCJ 



In the limit of small intergenerational mixing, P^vMsw ~~ * PivMSW i n such a way that Eq.i 
is well approximated by 

P*{E) w ^P 2 ,msw, (95) 

where P2 V msw is evaluated for Ara^ and 77 in a medium of effective density 

A e// « A cc + \a nc , (96) 

as inferred from Table [Til], with a level-crossing probability P# H governed by Ti e ff in Eq.(|36|). 
Given the respective standard v e <-> T and v e <-» v s effective densities: 

A ef f(i/ e «-> Vft tT ) = A cc , 
A e ff{Ve <-> v a ) = A C c + A NC , (97) 

Eq. ( |96"D immediately puts the adiabatic edge of the MSW transition between that of the two 
standard cases for the same oscillation parameters as shown in Fig. 01. 



On the non-adiabatic aspect of P2uMsw, we learn from Table [TT| that the mixing angle 
77 and the off-diagonal term Ai n d (both of which appear in the level-crossing probability) 
do not distinguish between the exact identities of vu and or equivalently, the sign of 
Am 2 +1 _. (They do, however, depend on the identity of v^t — a matter of interchanging 9 
and 0.) This is because the enhanced apparent mixing mode at Rh is that of V\ x «-> v^t- 
The state z/ lx , in turn, is a fixed, Am 2 +1 _ independent admixture of v\+ and i>i_, such that 
Anij +1 _ does not enter the scene so long as Eq. ([77]) is satisfied. The matter- induced mixing 
term varies with 9 and 0. Taking the case of v\ x <-> V2-, we observe that Ai n d vanishes 
exactly when 9 = 0. For 9 = 0, we have Ai n d = — ^p- sin 2rj = ' Al ^ c ' sin 2rj so that the ratio 
of matter-induced to vacuum mixing evaluated at resonance is given by 

A = = J0 3 2', « 0.04 -+ 0.14. (98) 

The contribution of matter-induced mixing (which is additive here by definition) is therefore 
relatively small. Given an effective density in Eq.(|96Q that is almost independent of 9 
and (provided they are small), approximately the same Am^ will provide a fit to the 
experimental data for any combination of 9 and 0. It then follows from Eq. fl98D that for the 
9 = case, sin 2 2rj lies in the range 

sin 2 2r] « (0.8 -> 0.9) sin 2 2cp, (99) 
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where sin 2tp is the fitted mixing parameter for the case where no matter-induced mixing is 



present (i.e., when 6 = 0). At the other extreme, Table III shows that for V\ x <-> v 2 -, matter- 
induced mixing takes full control if = 0. This mixing, however, is negligible. Following 
from Eqs.(|65|) to (|7|)> the quantity k 2 (x) that appears in the level-crossing probability is 
now proportional to A 2 nd , i.e., 

(| Awe I cos 2?? cos a sin #) 2 = (Am^) 2 sin 2 26 




cos 2 6 + 1 

« (8 x 10~ 4 -> 10" 2 )(Am 2 1 ) 2 sin 2 20. (100) 

Thus, given that sin 2 26 is sufficiently small, the level-crossing probability for (f> = is 
almost one. [If the resonance is sufficiently close to the centre of the sun where ^ ~ 2, 
matter-induced mixing may show itself by contributing to the level-crossing probability an 
equivalent of sin 2 26 x 10~ 2 by Eq.fllOOD.I 

We may carry out the same analysis for v\ x <-> v 2 +, which involves little more than 
interchanging 6 and <fi. 



C. Case C 

Case C comprises the following parameters: 

Am 2 +2 _ > Am 2 21 > Am 2 +1 _, Am 2 +1 _ ~ 10~ u eF 2 . (101) 

The mathematics that describes the resonant conversion of u e in the interior of the sun 
for this case is identical to that for Case B2 where the resonance Rh is one of enhanced 
v ix v 2 t mixing. Indeed, the two systems are physically identical. However, the suppression 
of v e <-> v' e oscillations is no longer solely a matter effect but is extended to the vacuum, 
leading to a vastly different phenomenology. Vacuum v e <-> v' e oscillations do not happen 
because of the extremely small Am 2 +1 _ which in turn corresponds to an oscillation length 
much longer than an astronomical unit. We may therefore treat the problem as though the 
lower resonance Ri (responsible for u e «-> u' e ) in Fig. [3] is completely absent. The v e survival 
probability is then given by Eq.(^0|) with Pz v msw replaced with P2uMsw, i>e., 

Pg(E) = \U e2h \ 4 + (1 - \U e2h \ 2 fP 2vMSW . (102) 

The MSW transition probability P 2 vMsw is that in Eq . (^2|) evaluated for the effective Hamil- 
tonian of Case B2 in Eq.Q85D with a squared mass difference Am 2 j^ = Am 2l and mixing 
angle rj defined in Eqs. (|87D and (|82|) respectively. Assuming small intergenerational mixing, 
Eq.( |102j ) reduces to 

P®(E) w P 2 „msw, (103) 

where P 2v msw contains non-standard features as described earlier for Case B2. Equation 
( |103| ) is represented graphically in Fig. 0. 
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IV. SOLUTIONS TO THE SOLAR NEUTRINO PROBLEM 



In this section, we locate the regions of parameter space that will give rise to the observed 
solar neutrino depletion for Cases Bl, B2 and C in an approximate way. Currently available 
experimental data suggest a significant depletion of the mid-energy neutrinos, hinting at a 
pre-defined shape for the v e survival probability. By going to regions of parameter space in 
which matter-induced mixing is small, the approximate v e survival probabilities in Eqs.(j54T), 
(B5|) and (1103) for Cases Bl, B2 and C respectively are very simply related to the standard 
two-flavour P2vMSW- Thus, by comparison with the standard solutions, we may gain a rough 
feeling for the necessary oscillations parameters for each case without performing an ab initio 
fit to the experimental data. 

Flux-independent data from Kamiokande and SuperKamiokande such as spectral distor- 
tion and day-night asymmetry provide yet another means to identify the allowed oscillation 
parameters |33H35|. Although the day-night effect is beyond the scope of this paper, we will 
be able to comment on the expected spectral distortion in Cases Bl, B2 and C. 



A. Case Bl 

Let us reiterate that the oscillation parameters for Case Bl are constrained by Eqs.( 
and fl45|) . For comparison, it is useful to define the ratio of the event rate with oscillations 
to that with no oscillations as 

n _ lE Pl(E)^(e,E)a(E)dE 

S^{e,E)a(E)dE ' ^ 

where (j)®(e,E) is the no-oscillation v e flux, a(E) is the detection cross-section, E is the 
experimental energy threshold, and P ee (E) is the v e survival probability averaged over pro- 
duction positions. By Eq . (|64|) , the energy-dependence of the v e survival probability for this 
case is contained entirely in the term P2 V msw- I n the small intergenerational mixing limit, 
we further approximate the adiabatic edge of P2 U msw as a step function 9(E^ — E), where 
E\ v is Ea defined in Eq.(^), such that 

P 2 umswU, = 9{Ef -E) + P%. (105) 

For comparison purposes, we also write down an analogous expression for a 2v system: 

P2»MSW\2» = &{E\ V — E) + Pg, (106) 



where 



_ 9 „ Ami cos2i7 

E\ v = 21 (107) 



2y/2G F N e 



Equation ( [LOT] ) comes from the resonance condition, A e ff = Acc = Am2!Cos2?7, for a 
standard v e v^ T system, and a comparison with Eq. (|63|) immediately leads us to 

2E\ V = E\ v . (108) 
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The scale heights 



1 dA eff 



that appear in the level-crossing probabilities Pr and Pr 
re almost exponential solar density profile plj. Thus, if we 



I A eff dt 

do not differ much due to t 
ignore matter-induced mixing by setting 9 = according to Table |I|, we may make the 
approximation: 

Pr = Pr, (109) 
such that, by Eqs.(H), (OH ~ B) and (TO) . 



Q iu (Am 2 , sin 2 2r)) ~ - + ^tt 2v (2Am 2 , ]- sin 2 2rf). (f fO) 

Here, fi^Am 2 , sin 2 277) denotes the value of Q in our 4z/ scheme evaluated for Am 2 and 
sin 2 2rj, and f22v(2Am 2 , | sin 2 2rj) is the ratio fl in the standard z/ e <-> z/ MT scheme evaluated 
at twice the squared mass difference and half the mixing. The rescaling of Am 2 and sin 2 2rj 
in the latter is dictated by Eq.( |109| ), such that the parameter Am 2 sin 2 2rj that is fed into Pj^ 



and P 2v respectively agree. Thus, in the extreme case of Am 2 +1 _ ~ 10~ 3 ey 2 , a minimum 
of \ of the original neutrino flux must be detected, while the maximum detectable flux is 
I as a direct consequence of maximal vacuum v e <->• v' e oscillations. The latter corresponds 
to the absence of MSW transitions, or equivalently, to Am 2 x and rj residing in regions of 
parameter space outside of that quoted in Eqs. (|50"D and ([H]) respectively. 

Equation (|110| ) allows us to virtually fit the experimental data by the use of existing the- 
oretical v e <->• z/ Mj t predictions for the various experiments in which smearing over production 
point and other energy dependence are already accounted for. Given a set of standard 2v 
MSW contours, we can pick out the necessary oscillation parameters for this case graphically 
by identifying each f2 2 ^ contour with a by Eq.( |110| ), and adjusting the Am 2 and sin 2 2rj 
scales accordingly. The area enclosed by the dot-dash line in Fig. ^ represents the region of 
parameter space in which the 2a bands (including both experimental and theoretical errors) 
of all five experiments overlap ||36||. Numerically, this region is defined by 



2 x 1CT 6 ~ Am 2 21 /eV 2 ~ 5 x fCT 5 , (111) 

and 

8 x lCT 4 ~ sin 2 2r/ ~ 0.1. (112) 



[Beyond sin 2 27/ ~ 0.1, Eqs. (|64|) and ( |110| ) become invalid as intergenerational vacuum oscil- 



lations increase in amplitude.] Note that our fitting procedure is approximate, but should 
nevertheless yield a solid indication of the 2a allowed parameter space. We also find that 
the la bands do not overlap in this case, suggesting that the solution is not acceptable at 
below 68% C. L.. Figure |7| shows the v e survival probability for several representative sets 
of oscillation parameters within the allowed region for Case Bl. 

A qualitative discussion of this virtual fit follows. Firstly, from Eq. (|17D , to Table |, large 
uncertainties in the Boron flux means that Kamiokande and SuperKamiokande's 2a bands 
necessarily span the entire region of MSW parameter space for this case. Secondly, the 
central values of the GALLEX and SAGE results are somewhat higher than the maximum 
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detectable ratio of | as predicted by the EPM. Hence, the lower bound on Am^ is set by 
the Gallium experiments, corresponding to the maximum amount of low- and mid-energy 
neutrinos that can be killed resonantly within limits. 

The central value of the Homestake result is ~ \ relative to SSM prediction, while at plus 
2cr, the ratio of measured to no-oscillation event rates does not quite reach the maximum 
of |. Thus, the upper bound on the squared mass difference Am^ and the lower bound on 
Amji sin 2 2r\ are determined by Homestake, representing the possible suppression patterns 
that the Boron spectrum may receive. 

Note that we have arrived at these assuming 9 = <p. A rough indication of the necessary 
regions of parameter space for cases of unequal 9 and can be obtained based on the analyses 
in Sec. [Hi]. In short, unequal 9 and will give more breadth to the allowed region. 



B. Case B2 

In Case B2, the parameters are constrained by Eq. flT7|) . The v e survival probability for 



this case in Eq . (|92|) consists of negligible constant terms due to averaged vacuum oscillations. 
This is a direct consequence of the strong suppression of v e <-> v' e oscillations in matter such 
that almost all of the original ^ e 's produced at p ^> pn will pass through a resonance. In 
this case, P®(E) plunges down to as low as ~ immediately to the high-energy end of the 
adiabatic edge according to Eq. ( p5| ) and Fig. f|. The correspondingly deeper pit in PivMsw 
implies that the adiabatic edge needs to occur at an even higher energy in order to maximise 
the number of low- and mid-energy neutrinos to be detected. This is attainable by choosing 
an even higher squared mass difference. 

The naturally smaller-than-standard effective density in Eq.Q^) places the adiabatic 
edge of Case B2 at a somewhat higher energy. The exact location of the edge relative to a 
standard edge, however, cannot be simply quantified since N e (x) and N n (x) do not exactly 
track each other in the sun. However, expressing A:// of Eq.([56|) in terms of Acc for the 
relevant solar densities, 

7 23 

4,// « ( g - ^Moo, (113) 

the shift of the adiabatic edge relative to the standard location is negligible. On the non- 
adiabatic side, if we neglect matter- induced mixing (by setting 9 = for v\ x «-> v^-, or 
= for v\ x «-> ^2+), Pr evaluated for this case will be approximately the same as P^ for 
a standard 2v system with the same Am 2 and 77. Hence, PivMsw\<±v ~ PzvMSwlzv an d 

Q Au (Am 2 , sin 2 2r)) w ^ 2 „(Am 2 , sin 2 2r/), (114) 



by Eqs.(|95D and (|104|) , and the symbols carry the same definitions as before. Thus, utilis- 
ing established 2v MSW transition probability contours, we obtain the allowed oscillation 
parameters for Case B2 by Eq.( |114j ) in a region defined by 



1(T 5 ~ Am 2 21 /eV 2 ~ 10" 4 , (115) 

and 
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1(T 4 ~ sin 2 2r/ ~ 4 x 1CT 3 , 



(116) 



including both experimental and theoretical errors at 2a as shown in Fig. ^|. This is a 
considerably smaller region than that for Case Bl. In particular, the Kamiokande and 
Super Kamiokande results now place an upper bound on the allowed sin 2 2rj such that we 
do not wipe out too many Boron neutrinos by the MSW mechanism in the non-adiabatic 
branch. The v e survival probability for this case for several representative sets of oscillation 
parameters within the allowed region is shown in Fig. [8]. 

C. Case C 

Case C is described by the parameters given in Eq. (|101|) . In view of Fig. ||, the similarity 
between the energy-dependences of this case and of the standard 2v (both v e <-> u^ T and 
^ e v s ) scenarios shows that due to the absence of vacuum v e <-> v' e oscillations, the 
necessary oscillation parameters for Case C will lie in between those of the two standard 
cases (see Table [IV]). The solution at 95% C. L. is shown in Fig. |6|. 

V. PROSPECTS FOR FUTURE EXPERIMENTS 

We now briefly discuss the implications of the various solutions on experimental observ- 
ables. In particular, we shall look at the Boron energy spectrum observed by the high-energy 
scattering experiments, the Beryllium line that will be detected by BOREXINO, and the 
charged to neutral current event rates to be measured by SNO. 

A. Boron spectral distortion 

Qualitatively, the underlying maximal v e <-> z/ vacuum oscillations in Cases Bl and B2 
will lead to recoil electron energy spectra for SuperKamiokande and SNO respectively that 
are almost flat with respect to SSM predictions. This can be seen by comparing the various 
survival probabilities in Figs. and || For Case Bl (Fig. |7|), the energy- independent \ that 
provides a lower limit to the flux depletion considerably softens the energy- dependence of 
the non-adiabatic branch of the v e survival probability. For Case B2 (Fig. §), the slope of 
the non-adiabatic branch is scaled down by a factor of two, relative to the standard MSW 
solutions, because of Eq.fl95|). For Case A, complete energy-dependence means that spectral 
distortion is absent as in the no-oscillation case. Precise predictions for the amount of 
deformation in energy-dependent cases relative to standard expectations cannot be obtained 
without performing an ab initio numerical fit, because of the vastly different oscillation 
parameters involved. However, it may be said for certain that the spectral distortions in 
Cases Bl and B2 are significantly weaker than those predicted by all other minimal v e <-> T 
and/or v e <-> v s schemes currently in the market ||38|| , and are somewhat stronger than the 
standard large mixing angle (LMA) scheme. 

Based on the deviation of the averaged measured electron kinetic energy from its 
standard value, SuperKamiokande's flux-independent data to date do not distinguish be- 
tween the standard 2v small mixing angle, the large mixing angle and the no-oscillation 
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solutions within la f35f . From Figs. [7] and |8|, we expect the quantities resulting from 



Cases Bl and B2 respectively to take on some intermediate values, compared with the stan- 
dard SMA and LMA solutions. In this respect, Cases Bl and B2 are consistent with spectral 
data to date. Case A is also acceptable. (A small reduction in the allowed regions for Cases 
Bl and B2 displayed in Fig. §]may result from a rigorous consideration of existing spectral 
distortion data. If this were to occur, its cause would be the encroachment of the adiabatic 
edge above the energy threshold for Super Kamiokande, which is presently ~ 6.5 MeV. An 
inspection of Figs. |7| and || however, reveals that the adiabatic edges occur at less than 
6.5 MeV in our admittedly approximate fits.) 

Observational effects associated with Case C are similar to those studied in Ref. 0. 
This case will give a recoil electron energy spectrum that is similar to that predicted by 
the standard 2z/ cases. Currently available flux-dependent and flux-independent data do 
not distinguish between this case and the standard v e <-> v^ T and v e <-> v s SMA scenarios. 
However, the SMA solution is preferred by SuperKamiokande over the LMA solution based 
on spectral data analyses |j35|. To this end, Case C looks promising. 



B. The Beryllium line 



With the exception of Case C which predicts an energy-dependence that is similar to the 
standard v e «-> v^ T and v e «-> v s scenarios, maximal vacuum v e <-> v' e oscillations lead to a 
Beryllium line that must be detected at ~ \ to ~ | the no-oscillation rate by BOREXINO. 
These deductions come from an inspection of Figs. [7| and || In particular, if Case B2 is 
valid, the Beryllium flux will be almost exactly halved, independent of the intergenerational 
oscillation parameters (provided they lie within the allowed region shown in Fig. |6|). In this 
respect, Cases Bl and B2 are clearly distinguishable from the standard 2v SMA schemes. 
Needless to say, Case A being energy-independent will exactly halve the Beryllium flux. 



C. Charged to neutral current rate 



In Sec. [TI], we obtained a range of the ratio of charged to neutral current event rate at 
SNO for each of our several cases, based on flux- dependent data from SuperKamiokande 
alone. Having identified in Sees. [I Ij and |I^ the allowed shapes of the various v e survival 
probabilities which are now constrained by five experiments, we may narrow these ranges, 
using expressions developed in Sec. 0. We will approximate the energy-averaged v e survival 
probability (-P®) as the P®(E) evaluated at an energy of lOMeV for our various cases. The 
ratios are shown in Table M. 

Considerable overlapping between the ranges for Cases Bl, B2 and C means that 
is perhaps not the best experimental observable for their disentanglement. However, by 
measuring alone, these cases are clearly distinguishable from the standard 2v scenarios, 
both pure active and pure sterile, and from Case A. 
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VI. IMPLICATIONS OF LSND 



The claimed observation of z/ M — > v e and — > v e by LSND suggests small angle mixing 
between neutrino states separated by a squared mass difference of ~ 0.1 — > 10 eV 2 f39|. This 



feature can be easily incorporated into Case A. For Cases Bl, B2 and C, and adhering to the 
standard neutrino mass hierarchy, the required oscillation length suggests indirect <-> v e 
oscillations through a sufficiently heavy u T and/or v' T . However, these are quite suppressed 
if consistency with all other neutrino experiments to date is sought. Reconciliation with the 
LSND result occurs at about the 3a level in a minute region of parameter space (see Refs. 
jTJ and |4(| for relevant discussions). 

For greater consistency with LSND, the standard mass hierarchy m\± < m 2 ± < m 3 ± 
must be altered by interchanging the second and third generation neutrinos (inverse mass 
hierarchy between v^/v' and v T /v' T ). In this scenario, the MSW partners of v e are v T and/or 
v' T . The v,j,/v' pair is now placed at ~ 0.1 — > 10 eV 2 above the v e jv' e pair on the squared mass 
spectrum. With an appropriate mixing angle, this will lead to direct «-> v e oscillations, 
thereby accounting for the LSND result. 

For both of these scenarios, relic neutrino asymmetry generation via ordinary-mirror 
neutrino oscillations may not be strong enough to ensure consistency with Big Bang Nu- 
cleosynthesis ||41||. The reason is that the LSND Am 2 tends to be smaller than the Am 2 



values favoured by this mechanism (see the third paper in Ref. 123). However, a detailed 
calculation would need to be performed to be sure of this. If this mechanism fails, then con- 
sistency with Big Bang Nucleosynthesis can be obtained by postulating that a sufficiently 
large neutrino asymmetry (~ 10~ 5 eV 2 ) is created at a high temperature scale by some phys- 
ical mechanism unassociated with ordinary-mirror neutrino oscillations (Ref. |32| discusses 
this type of scenario in more detail). 



VII. SUMMARY AND CONCLUSIONS 

In order to account for the significant depletion of solar neutrinos measured by five ex- 
periments to date, we have invoked MSW-enhanced intergenerational mixing in addition to 
maximal vacuum v e <-> v' e oscillations prescribed by the EPM. Approximate analytical ex- 
pressions for the u e survival probabilities as functions of neutrino energy have been obtained 
for several possible neutrino mass hierarchies assuming small vacuum intergenerational mix- 
ing. These expressions were then compared with well-established 2u solutions to identify the 
approximate regions of parameter space that can simultaneously explain the apparent solar 
and atmospheric neutrino anomalies. (Note that only those parameter space regions which 
feature small matter-induced mixing were examined in depth. The approximate allowed re- 
gions plotted in Fig. [| assume this restriction. Some indication of the effect of non-negligible 
matter- induced mixing is discussed in Sec. |TTT| . ) The v e survival probabilities for the Cases 
A, Bl, B2 and C considered herein exhibit considerable differences from each other and from 
the standard v e <-> v^ T and v e <-> v s small and large mixing cases. 

These new, distinguishing features are, in principle, observable by future solar neutrino 
experiments, and have been briefly discussed. The results are summarised in Table V. 
By utilising Boron neutrino spectral distortion, the Beryllium neutrino flux and the ratio of 
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charged to neutral current event rates, all the EPM possibilties can be distinguished from the 
standard 2v MSW solutions, and, with the exception of Case C, from models incorporating 
only one sterile neutrino. The four cases within the EPM also yield different outcomes, 
except for Cases Bl and B2 which exhibit an overlap for these solar neutrino observables. 
Fortunately, Cases Bl and B2 can be differentiated through atmospheric neutrino data. 
The large Amf +l _ that defines Case Bl leads to significant v e <-> v' e atmospheric neutrino 
oscillations. This case is actually already disfavoured by the recent atmospheric neutrino 
data from SuperKamiokande (see Ref . [[13[ for detailed discussions) . Note that Case C within 
in the EPM is very similar to the scenario analysed in Ref. 0. 

We should also note that the implications of MSW solutions within the EPM for the day- 
night effect have yet to be examined. In addition, one relevant region of parameter space 
has not been explored in detail in this paper: Am 2 +1 _ values in the approximate range 
10~ 6 — > 1CT 4 eV 2 which are intermediate between Cases Bl and B2. In this regime, the 
mass-squared difference between v 1+ and v\- is comparable to the MSW intergenerational 
mass difference, and our approximation scheme is no longer reliable. This region is perhaps 
best explored numerically, a task beyond the scope of this paper. 

The Exact Parity Model is, in part, an explicit theory of light, effectively sterile, neu- 
trinos. Its characteristic ordinary-mirror neutrino maximal mixing feature receives strong 
experimental support from the atmospheric neutrino data. Various possibilities for solv- 
ing the solar neutrino problem by either averaged vacuum v e z/ oscillations alone (Case 
A), or several amalgams of MSW-enhanced and vacuum oscillations (Cases Bl, B2 and C), 
exist within the EPM. Future solar neutrino experiments should narrow the possibilities 
considerably. 
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TABLES 



TABLE I. Solar neutrino measurements and theoretical expectations within the standard so- 
lar model of Bahcall and Pinsonneault (SSM-BP) (1995) pi] . Capture rates for Homestake [14], 



GALLEX U and SAGE QlO|] are given in SNU, where 1 SNU = 10~ 36 capture per atom per second. 
For Kamiokande [12] and SuperKamiokande O], the measured neutrino flux is given in 



10 6 cm 2 s . The associated statistical and systematic errors (lex) are quoted for each experiment. 

Experiment Measurement SSM-BP 

Homestake 2.56 ± 0.16 ± 0.14 9.31^ 

GALLEX 69.7 ± 6.7±1:§ 1371? 

SAGE 72±ig±| 1371? 

Kamiokande 2.80 ± 0.19 ± 0.33 6.62±5;^ 

SuperKamiokande 2 -5llai3 ± 0.18 6.62±2;^ 



TABLE II. The mixing angle rj, the effective density A e ff and the off-diagonal matter-induced 
mixing term Ai n d evaluated for various combinations of v\h and u^i for Case Bl. 

v\h <-> V2l sin?? Aeff And 

vi- vi- sin (f> 

v\- <-> V2+ —j — = — 5= ^(cos^ + sin 2 ^) + A NC sm2r]sm{e - (f>) A NC cos 2rj sm(6 - (f>) 

y cos 2 </>+sin 8 

v\+ «-> V2+ sin 9 ^Np 

v\+^v 2 - i sin0 g (cos 2 6 + sin 2 0) - A JV csin2r/sin(6' - 0) -A^c cos 2r?sin(6» - 0) 

Y cos 2 0+sin </> 



TABLE III. The mixing angles ip and rj, the effective density A e ff and the off-diagonal mat- 
ter-induced mixing term Ai n d evaluated for various combinations of uu, v\% and v^t for Case B2. 

U ^v 2i sin-0 sin?? A eff A ind 
v\h 

l (l + cos 2 



■ ■ ^ / «** t +^^9^9,/,™^-^ -^c[cos2r / cos^ S in( 



y^os 2 0+cos 2 </> Vl+cos 2 e 



-2 sin 2r] cos sm(6/ — <p)\ z 



^2- 



— 9^(1 + COS 2 i/i . ... . (/] 

+^[cos2r7sin2^cos^-^ -^w[cos 277sm^sm(0 - 4>, 

yWe+cos^ v^I+c^ + 2 [cos2r/sm2^cos^ 0j + l sin2 in2 ^, cos( 0_ 0)] 

— 2sin2?7sin^sin(6' — <p)\ 

AqC (1 + COS 2 (/)) 

vi ~ - - ■** -r^_ + V [«* ^ Bin 2^ cos(e - ^) A r [co s 2,7 co ; ^ ~j3 

z+ vW e+cos 2 4> yT+W^ 2 l / y v v; -i s i n 2r7sin2?/)cos(0- 0) 

1 2 sin 2r7 cos ?/> sm(# 11 - 



^(i + cos 2 ^) " : . . , . ~ ~r 

« . n , fn ,\ Anc cos 277 sm to sin (0 — d>) 

*' ^^2+^ — = — cos27/sm2V>cos(0- <£) V ' , , v „ 

V\- Ycos 2 e+cos 2 </> yi+cos 2 </) „._.,.//, ,m — n sm 2t7 sin 2w cos(# — <p)\ 

+2sm 2r] smip sm(0 — (p)\ z 



cos 4> sin 
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TABLE IV. The allowed intergenerational oscillation parameters Am 2 and sin 2 2r\ for Cases 
Bl, B2 and C. The best fit oscillation parameters in the small mixing angle (SMA) and the large 



mixing angle (LMA) solutions for standard 2v scenarios [33,371] are also included. 



Scheme 



Am 2 /eV 2 



5 x 10" 6 
1.6 x 1(T 5 

4 x io~ 6 

2 x 1(T 6 

i(r 5 

4 x 1(T 6 



sin 2 2rj 



8 x 1(T 3 
0.63 
10~ 2 
8 x 10" 4 -> 0.1 
10~ 4 -> 4 x 10~ 3 
8 x 10~ 3 10~ 2 



Standard v e 



SMA 
LMA 



Standard v e *-> f MjT 
Standard u e <-> z/ s SMA 
Case Bl 
Case B2 
Case C 



5 x 10" 5 
10" 4 
5 x 10" 6 



TABLE V. Predictions for future experiments. The amount of Boron spectral distortion as- 
sociated with each case is qualitatively compared with the standard 2v predictions (maximal for 
SMA and minimal for LMA). The approximate Beryllium fluxes to be measured by BOREXINO 
prescribed by our various cases relative to the no-oscillation flux are also compared. The last 
column shows the predicted ranges of the ratio of charged to neutral current event rate relative to 
the no-oscillation rate at SNO. 

Scheme Boron spectral distortion Be flux ( pr gg^ cd ) r d for SNO 

Standard v e *-* Va T SMA maximal ~ 0.25 — > 0.4 

Standard v e *-* v„ T LMA minimal ~ 0.4 ~ 0.2 

Standard u e <-> v s SMA maximal ~ 1 

Case A none \ 1 

Case Bl intermediate \^\ ^.5 —>■ 0.75 

Case B2 intermediate ~ \ 0.4 — >■ 0.7 

Case C maximal ~ 0.45 — > 0.6 
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FIGURES 

FIG. 1. Level-crossing diagram for Case Bl. The labels v™, where i = 1£, lh, 21, 2h, denote 
the instantaneous mass eigenstates. The letter R labels the intergenerational MSW resonance. 

FIG. 2. The v e survival probability at the Earth for Case Bl with Am|i = 5 x W~ 6 eV 2 and 
sin 2 2rj = 8 x 1CP 3 (solid line) for v e produced at the centre of the sun. The dashed and dotted lines 
represent, respectively, the survival probabilities for the standard v e <-> i/ MiT and v e <-> v s scenarios 
evaluated for the same oscillation parameters. 

FIG. 3. Level-crossing diagram for Case B2. The labels u™, where i = 1£, lh, 2£, 2h, denote 
the instantaneous mass eigenstates. Two resonances are identified and labelled as Rl and Rh, 
where Rh is the resonance of interest. 

FIG. 4. The v e survival probability at the Earth for Case B2 with Am^ = 5 x W~ G eV 2 and 
sin 2 2r\ = 8 x 10~ 3 (solid line) for u e produced at the centre of the sun. The dashed and dotted lines 
represent, respectively, the survival probabilities for the standard v e *-* f MiT and v e *-* v s scenarios 
evaluated for the same oscillation parameters. 

FIG. 5. The v e survival probability at the Earth for Case C with Amfi = 5 x l^T^eV 2 and 
sin 2 2r\ = 8 x 1CP 3 (solid line) for u e produced at the centre of the sun. The dashed and dotted lines 
represent, respectively, the survival probabilities for the standard u e *-* Va tT an d v e v s scenarios 
evaluated for the same oscillation parameters. 

FIG. 6. The approximate allowed regions for Cases Bl (area enclosed by the dot-dash line) 
and B2 (area enclosed by the dashed line). These correspond to the regions in which the 2a bands 
(including theoretical uncertainties) of all five solar neutrino experiments (see Table |J) overlap, 
respectively for Cases Bl and B2. Note that these regions will differ slightly from those generated 
from a x 2 -analysis. The allowed region for the standard 2v case at 95% C.L. (dotted line) is shown 



here for the purpose of comparison [33]. The solution to Case C is similar to this. 



FIG. 7. The v e survival probability at the Earth for Case Bl evaluated for various Am|i and 
sin 2 2r] shown on the graphs (solid line) for v e produced at the centre of the sun. These parameters 
lie within the allowed region for Case Bl but are not necessarily the best fit parameters. For the 
purpose of comparison, the survival probabilities for the standard v e *-* f„ jT small angle solution 



(dashed line) and large angle solution (dotted line) (see Table IV for the best fit parameters) are 
also shown. 

FIG. 8. The v e survival probability at the Earth for Case B2 evaluated for various Am^i and 
sin 2 2r] shown on the graphs (solid line) for v e produced at the centre of the sun. These parameters 
lie within the allowed region for Case B2 but are not necessarily the best fit parameters. For the 
purpose of comparison, the survival probabilities for the standard v e *-* z/„ jT small angle solution 



(dashed line) and large angle solution (dotted line) (see Table |V| for the best fit parameters) are 
also shown. 
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